Spontaneous Emergence of Persistent Spin Helix from Homogeneous Spin Polarization 
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We demonstrate that a homogeneous spin polarization in one-dimensional structures of finite 
length in the presence of Bychkov-Rashba spin-orbit coupling decays spontaneously toward a persis- 
tent spin helix. The analysis of formation of spin helical state is presented within a novel approach 
based on a mapping of spin drift-diffusion equations into a heat equation for a complex field. Such 
a strikingly different and simple method allows generating robust spin structures whose properties 
can be tuned by the strength of the spin orbit interaction and/or structure's length. We generalize 
our results for two-dimensional case predicting formation of persistent spin helix in two-dimensional 
channels from homogeneous spin polarization. 
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The helical wave of rotating spin orientation is re- 
ferred to as the spin helix. There is a significant interest 
to spin helix configurations in semiconductor materials 
since the electron spin relaxation of such spin configura- 
tions can be partially [1-3 J or even completely suppressed 
[4j |5] . While a partial suppression of spin relaxation in 
two-dimensional systems becomes possible in the pres- 
ence of only Bychkov-Rashba |6J spin-orbit coupling (see 
Refs. [TJ [2J) ; the complete suppression of spin relaxation 
requires a specific combination of Bychkov-Rashba and 
Dresselhaus |7J interactions as it was demonstrated in 
Ref. [4]. More generally, the relaxation of the spin he- 
lix is an example of situations [THSJ I5HT7] when electron 
spin relaxation scenario deviates from the predictions of 
D'yakonov-Perel' theory |18j . 

Experimentally, the spin grating technique [19J is typ- 
ically used [3j |5] to create spin helical configurations in 
semiconductors. In this method, a sample is illuminated 
by a pair of pump beams with orthogonal linear polariza- 
tions. The interference of such beams results in a spacial 
modulation of light helicity. Correspondingly, through 
the optical orientation effect, a modulation of spin polar- 
ization in the form of spin helix is produced. Moreover, a 
spin injection from a ferromagnetic material into a semi- 
conductor can also be used to excite a spin helix [20J. 
In this approach, the rotating spin polarization is caused 
by coherent spin precession of electrons drifting in an 
applied electric field. However, the present authors are 
not aware about any experimental studies of spin helixes 
excited by spin injection. 

In this Letter, we propose an alternative approach 
to induce spin helical configurations. Specifically, we 
demonstrate that in one-dimensional (ID) systems of fi- 
nite length with Bychkov-Rashba spin-orbit coupling the 
spin helical configurations emerge in the process of re- 
laxation of homogeneous spin polarization (see Fig. [TJ. 
Mathematically, such a strikingly unexpected transfor- 



mation of homogeneous spin polarization into the per- 
sistent spin helix occurs when we introduce boundary 
conditions on electron space motion to describe finite- 
length structures (in infinite systems the homogeneous 
spin polarization decays exponentially as predicted by 
D'yakonov-Perel' theory |18|). Using a novel approach 
that maps spin drift-diffusion equations into a heat trans- 
fer equation for a complex field we find the exact time- 
dependence of the spin polarization dynamics. It is inter- 
esting that the amplitude of the resulting spin helix has 
an oscillatory dependence on the system's length. Below, 
we provide an intuitive explanation of this result based 
on properties of solution of heat equation. Moreover, it 
is necessary to emphasize that our theory is generalized 
for the case of two-dimensional (2D) channels and can 
be straightforwardly verified experimentally. In partic- 
ular, experimentally, the homogeneous spin polarization 
can be easily created using the optical orientation by cir- 
cularly polarized light. Therefore, we believe that our 
approach would simplify tremendously the generation of 
long-living spin helical configurations in semiconductor 
structures and advance the field of spin storage in semi- 
conductors. 
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FIG. 1: (Color online) Schematics of spontaneous transfor- 
mation of homogeneous spin polarization into persistent spin 
helix in a finite length system with Bychkov-Rashba spin-orbit 
coupling. 
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FIG. 2: (Color online) Dynamics of formation of persistent spin helix from homogeneous spin polarization pointing in z 
direction at t = 0. These plots were obtained using Eq. ( |13[ ) at r\L — 15.45. This value of the parameter r\L corresponds to 
the second local maximum of spin helix amplitude shown in Fig. [3] 



Let us consider dynamics of electron spin polariza- 
tion in a ID system of a length L in x direction in the 
presence of Bychkov-Rashba spin-orbit coupling. In one- 
dimensional limit, spin drift-diffusion equations [2 J can 
be written as 
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where D = £ 2 /r is the coefficient of diffusion, A = 
d 2 /dx 2 , C = 2r]D is the constant describing spin rota- 
tions, 7 = r] 2 D/2 is the coefficient describing spin relax- 
ation, r] = 2amh~ 1 is the spin precession angle per unit 
length, a is the spin-orbit coupling constant, m is the ef- 
fective electron mass, £ is the mean free path and r is the 
momentum relaxation time. It follows from Eq. Q that 
y component of spin polarization, S y , is not coupled to 
any other component of spin polarization. Consequently, 
selecting S y (x, t = 0) = we can safely take out S y from 
our consideration. Eqs. ( T|3 > are complimented by stan- 
dard boundary conditions [21 J 
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Here, Y = [x = 0,x = L). Mathematically, the bound- 
ary conditions Q are so-called third-type boundary con- 
ditions. This specific form of boundary conditions con- 
serves the spin polarization of electrons that scatter from 
the sample edges. We assume that at the initial moment 



of time the spin polarization is homogeneous and points 
in z direction, that is 



S x (x : t = 0) = 0, S z (x,t = 0) = S - 



(5) 



Let us introduce a complex polarization S = S x + iS z . 
It is straightforward to show that Eqs. ( T|3 ) and bound- 
ary conditions Q can be rewritten in a more compact 
form using S: 
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Defining a complex field u(x,t) by the relation 



u(x, t) 



(8) 



we find that Eq. (|6| transforms into the heat equation 
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dt dx 2 ' 

supplemented by Neumann (or second-type) boundary 
conditions 
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Moreover, it is worth noticing that the initial conditions 
for u(x,t) are related to the initial conditions for S as 



u(x,t = 0) =e~ ir]X S(x,t = 0). 



(11) 



Consequently, the initially homogeneous spin polariza- 
tion in z direction (Eq. (|5|) corresponds to a spatially 
modulated complex field 

u(x, t = 0) = Sq sin (rjx) + iSo cos (rjx) . (12) 
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The solution of Eq. Q with the boundary condi- 



tions (10) and initial condition (12) was obtained by the 
method of separation of variables. It can be presented in 
the form 



S(x,i) _ .sin(r/L/2) 
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This is our main analytical result describing dynamics 
of spin polarization in ID finite-length structures. Note 
that S x and S z components of spin polarization are given 
by real and imaginary parts of Eq. (13), respectively. 



The first term in the right-hand side of Eq. (13) de- 



scribes the persistent profile of spin polarization (in the 
form of spin helix) emerging at long times. Concerning 
the second term in the right-hand side of Eq. ( fl3| ), it 
governs the dynamics of transformation of the initially 
homogeneous spin polarization into the persistent spin 
helix. Fig. [2] demonstrates dynamics of S z and S x com- 
ponents of spin polarization given by Eq. (13). It is 



clearly seen that the initially homogeneous spin polariza- 
tion in z directions transforms into the persistent spin 
helix with an (infinitely) long lifetime. 

Explicitly, in the long time limit, the spin polarization 
is given by 

S x (x, t = +oo) = -g Sm( ^{ 2) sm(rj(x - L/2)), (14) 
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In these equations the factor sm(r)L/2)/(rjL/2) defines 
reduction of the spin helix amplitude with the respect to 
the initial amplitude of homogeneous spin polarization 
So. We plot this function in Fig. [3] It is interesting 
that the spin helix amplitude is an oscillating function of 
the parameter rjL and takes zero values when rjL = 2ixn 
where n is a positive integer. The positions of local max- 
ima can be found numerically. In particular, positions 
of four local maxima shown in Fig. [3] are 8.987, 15.450, 
21.808, 28.132. 

The heat equation is the best starting point to under- 
stand the oscillatory dependence of spin helix amplitude 
on rjL depicted in Fig. [3j Accordingly to Eq. ( [TT] ), the 
initially homogeneous initial condition (Eq. ([5| for spin 
diffusion equations transforms into a modulated initial 
condition for the heat equation. As the solution of heat 
equation in the given context represents simply the pro- 
cess of temperature equilibration along the system, an 
integer number of modulation periods results in zero av- 
erage "temperature" and, correspondingly in zero spin 
helix amplitude. Moreover, we would like to mention 
that the spin helix formation process is described by a 
series of exponentially decaying terms whose time con- 
stants are given by r n = L 2 / (tt 2 n 2 D) . The longest of 
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FIG. 3: (Color online) Normalized amplitude of the persistent 
spin helix as a function of r\L. Insets show schematically 
distributions of S x and S z at several specific values of rjL as 
indicated by arrows. Positions of minima and maxima points 
of the amplitude are discussed in the text. 



these times t\ = L 2 /(tt 2 D) provides the time scale of 
the transformation process. The dependence of spin pro- 
cess on L is intuitively clear as electrons should "feel" 
the system's length before the transformation ends. It's 
also interesting that such a time can be longer or shorter 
then the relaxation time of homogeneous spin polariza- 
tion Th = l/(Dr] 2 ). In particular, n/r^ = (tjL/tt) 2 mean- 
ing that n < Th when r/L < the times are the same 
when rjL = 7r, and t\ > when rjL > tt (see also Fig. 

In order to obtain an additional insight on spin re- 
laxation of the radial spin helix, we have performed ex- 
tensive Monte Carlo simulations employing an approach 
described in Refs. |8J and [22\. This Monte Carlo sim- 
ulation method uses a semiclassical description of elec- 
tron space motion and quantum-mechanical description 
of spin dynamics (the later is based on the Bychkov- 
Rashba coupling term). All specific details of the Monte 
Carlo simulations program can be found in the references 
cited above and will not be repeated here. A spin con- 
servation condition was used for electrons scattering from 
system boundaries. Generally, all obtained Monte Carlo 
simulation results are in perfect quantitative agreement 
with our analytical predictions thus confirming the mech- 
anism of formation of persistent spin helix from homoge- 
neous spin polarization. A comparison of selected ana- 
lytical and numerical curves is given in Fig. [4j 

The results reported in this paper can be readily gen- 
eralized for the persistent spin helix in two dimensions 
[4]. Indeed, it can be easily seen that in the case of 
equal strength of Bychkov-Rashba and Dresselhaus spin- 
orbit interactions, a = j3 (where j3 is the Dresselhaus 
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FIG. 4: (Color online) Long-time distribution of S z at r/L = 
8.987 found employing Monte Carlo simulation approach. 
The analytical curve is obtained using Eq. (15). The Monte 
Carlo simulation was performed for 10 5 electrons in GaAs 
structure of 1.7/im length. This plot obtained using the pa- 
rameter values r = O.lps, / = lOnm, a = 3-10 _12 eV m. Inset: 
orientation of 2D channel for two-dimensional spin helix ex- 
citation experiments. 



spin-orbit coupling constant), the equations of spin dif- 
fusion in 2D [4J take the general form of Eqs. ([l][3]). 
Therefore, introducing appropriate boundary conditions, 
namely, reducing the system into a 2D channel in [-110] 
direction (see the inset in Fig. |4|, we obtain the situa- 
tion completely equivalent to that in ID from the point 
of view of spin dynamics. Taking into account recent ex- 
perimental demonstration of persistent spin helix |5J the 
emergence of persistent spin helix from homogeneous spin 
polarization can be straightforwardly detected. Finally, 
we would like to note that the amplitude of persistent 
spin helix can be increased by a repetitive excitation of 
homogeneous polarization by a train of laser pulses. 

In summary, we have demonstrated that persistent 
spin helix forms in the process of relaxation of homo- 
geneous spin polarization in finite length systems. This 
observation can be used as a different technique for creat- 
ing spin helical structures in semiconductors. The solu- 
tion of spin drift-diffusion equations describing formation 
of persistent helix was derived analytically and numeri- 
cally using Monte Carlo simulation approach. The re- 
sults obtained in both ways are in perfect agreement. It 
is interesting that the persistent helix amplitude demon- 
strates an oscillatory dependence on the system length 
and strength of spin orbit interaction. Therefore, the con- 



trol of spin helix characteristics is achievable via appro- 
priate choice of the above mentioned parameters. This 
suggested technique facilitates generation of spin helical 
states and can be used in both one- and two-dimensional 
geometries. 
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